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To reduce errors in light-quark mass determinations, it is now necessary to consider electromag- 
netic contributions to light-meson masses. Calculations using staggered quarks and quenched 
photons are currently underway. Suitably-extended chiral perturbation theory is necessary to ex- 
trapolate the lattice data to the physical limit. Here we give (preliminary) results for light-meson 
masses using staggered chiral perturbation theory including electromagnetism, and discuss the 
extent to which quenched-photon simulations can improve quark-mass calculations. 
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1. Introduction 

As fundamental parameters of the Standard Model, quark masses are important to high-energy 
phenomenology. Lattice gauge theory can be used for precision, non-perturbative calculations of 
these masses, with state-of-the-art calculations using a vacuum including three or even four light, 
dynamical quarks: and c. In addition, some recent calculations have included electromag- 

netic effects via a quenched photon field. This has been done, for example, by the RBC Collabora- 
tion [1]. In calculations of quark masses by the MILC Collaboration [2], however, electromagnetic 
effects have in general been estimated based on continuum calculations. For the calculation in 
Ref. [2], the uncertainty due to electromagnetic effects dominates m M /m^ and is significant for m u 
and mj, as shown in Table 1. To put this error on better footing and to gain the higher precision 
desired by phenomenologists, the MILC collaboration is working to reduce the electromagnetic 
uncertainty in m u /md [3, 4]. The "rooted staggered chiral perturbation theory" (rSZPT) expres- 
sions reported in this work are a necessary step towards that goal. Our calculations do not assume 
quenched photons but allow the sea quark charges to take nonvanishing values. For comparison 
with simulations with quenched electromagnetism one can simply set the sea quark charges to 
zero. 

2. Electromagnetism and Quark Masses 

Quark masses are obtained by tuning bare-quark masses such that lattice-calculated meson 
masses match experiment. Taking the K + mass as an example, one tunes the QCD-only lattice 
result (M^ + )qcd to match experiment with electromagnetic effects subtracted off [2]; 

(m| + )qcd = M 2 K+ - (M 2 K± -M£ )em. (2.1) 

where is the experimental value, and the subscript "EM" denotes the purely electromagnetic 
contribution 1 . 

To obtain a value for (M^ ± — M^ )em, one uses an observation made by Dashen in the late 
1960's that the leading order electromagnetic effects are identical in the kaon and pion systems [5]. 
Specifically, 

(Mf ± -M 2 k0 )em = (M 2 K± -M 2 k0 ) em . (2.2) 

Corrections to Eq. (2.2) are referred to as the "violation of Dashen's theorem" and are often 
parametrized as 

(M| ± -M 2 k0 )em = (l+A E )(Ml± -M^ )em. (2.3) 

A lattice calculation of the ratio (M^± — M^ )em/(M 2 ± — M 2 )em yields a value of 1 + Ae in 
which common systematics cancel. This can be combined with experimental measurements of 
(M 2 + — M 2 )em to obtain the kaon splitting needed for Eq.( 2.1). 

In Ref. [2], A E was based on estimates from continuum calculations. A direct lattice calcula- 
tion of Ae will allow for improved accuracy and precision. Simulations by the MILC collaboration 
using "asqtad" staggered quarks [6] are reported in these proceedings [3]. To complete the calcu- 
lation, the rSZPT calculation for the meson mass must also be done. We report our preliminary 
results for that calculation [7] here. 

l We assume, for now, that the electromagnetic effects are small for MjL. 
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m u 


md 


m u /m d 


central value 


1.9 


4.6 


0.42 


statistics 


0.0 


0.0 


0.00 


lattice-systematic 


0.1 


0.2 


0.01 


perturbative 


0.1 


0.2 




electromagnetic 


0.1 


0.1 


0.04 



Table 1: Quark masses as given in Ref. [2]. Electromagnetic effects are the leading error in the ratio m M /mj 
and significant in m u and mj. Electromagnetic errors are estimated from continuum calculations. Details 
about this and the determination of other errors can be found in Ref. [2] 



2.1 The Neutral Pion 

The true 7T° is a uu — dd state. Because of electromagnetic disconnected-diagram contributions 
to its mass, even for m u = md, it is difficult to simulate. Instead, one simulates a meson created from 
a light quark and a light anti-quark that are distinct flavors and have opposite electric charges [1,3]. 
We call this neutral meson the 7i f , and make the replacement 

Ml±-Ml ^Ml±-Mi, (2.4) 



for the mass splitting. In the continuum, this approximation neglects some electromagnetic contri- 
butions to the physical 7T° 

Ho ~Ml,) EM = ^j^Ml (lnMl + l)+e 2 Ml (LECs) , (2.5) 

where "LECs" stands for a linear combination of the low energy constants. At the physical point, 
the difference given in Eq. (2.5) is small compared to the electromagnetic correction to the charged 
pion, <5 EM M^±, which is dominated by the kaon mass: 

8 EM M 2 n± * ^^MflnMl + ^MKLECs) , (2.6) 



with 

4e_ 2 C 
7 

where e is the fundamental electric charge, and C is a leading order (LO) LEC 



A EM = ^^, (2.7) 



3. Overview of the Calculation 



We split the overview into two parts. Section 3.1 reviews the needed results from previous 
calculations, while Sec. 3.2 describes new work. 
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3.1 Previous Calculations 

The leading-order, continuum, chiral Lagrangian including electromagnetic effects, but drop- 
ping the pure gauge-field terms, is [8] 

jgf(2) = l -f(d^d,L) + l -ftfz + X tf) - ±m 2 (<S>) 2 + e 2 C(QLQtf), (3.1) 

where we work in Minkowski space, for ease of comparison with Ref. [8]. Here, / « 130 MeV is 
the leading-order pion decay constant in the chiral limit, and % = 2Bo^, where M is the usual 
quark-mass spurion. 

The first two terms in Eq. (3.1) are the standard kinetic and mass terms. The meson fields are 
contained in 

U 7T+ K + \ 

£ = exp(i20//), = I 7T _ D K° 



K~ K° S 



(3.2) 



where diagonal entries of are comprised of u,d,s quark-anti-quark pairs; e.g. U = uu. 

The anomaly term ^m^O) 2 gives mass to the T]' . It also causes the "mesons" on the diagonal 
of <I> to have disconnected propagators [9, 10]. Usually, one takes mo — » 00 at the start, decoupling 
the r\' . For a partially quenched and/or staggered calculation, though, it is more convenient to keep 
the anomaly term and use the simple U, D, S basis along the diagonal of <&. In the end, we do take 
mo —> 00 [11]. 

Electromagnetic effects enter via the the photon field A M in the covariant derivative 2 

d^Z = d^Z-ieQA^Z + iZeQA^. (3.3) 

and the term e 2 C(QLQl^), where C is determined in terms of the pion splittings via Eq. (2.7). Q is 
the quark (electric) charge matrix 

£> = diag(q u ,q d ,q s ) (3.4) 

and has the property (Q) = 0. With dim(eQ) = dim(p), dim(A M ) = 1, and e ~ p, the overall 
power-counting scheme is p 2 ~ M 2 e 2 . 

For lattice calculations it is useful to consider the partially quenched case where valence and 
sea quarks are distinct, and may have different mass. Taking the limit where valence and sea quark 
masses are equal, known as full QCD, recovers usual chiral perturbation theory (ZPT) results. We 
refer the reader to the literature [9, 10, 12] for details. For studying electromagnetic effects, the 
partially-quenched structure allows one to separate valence and sea quark charge contributions in an 
advantageous way. The partially-quenched ZPT calculation including electromagnetism has been 
done by Bijnens and Danielsson [13]. This is extremely useful and provides us with cross-checks 
for the staggered calculation. 

The rSZPT calculation for meson masses done by Aubin and Bernard [14], without electro- 
magnetism, addresses the complexities confronted when dealing with staggered quarks. For each 
continuum quark simulated, there are four staggered quarks. These are distinguished by a quantum 
number referred to as "taste". When combined to form mesons, the result is a set of sixteen mesons 



2 Vector and axial vector current sources have been set to zero, since they are not needed here. 
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(a) (b) (c) (d) 

Figure 1: Feynman diagrams that contribute to the meson-mass at 0(p 4 ). Straight lines are the pseudoscalar 
meson propagator and wiggly lines are the photon. A filled dot represents an jSf ^ vertex, while an open 
square represents an ^ insertion, (a) photon tadpole; (b) photon sunset; (c) meson tadpole; (d) 0(p 4 ) 
tree-level insertion. 

of different tastes, which, at 0(a 2 ) in the lattice spacing a, can be grouped into five irreducible 
representations labeled: pseudoscalar, axial- vector, tensor, vector, singlet. The pseudoscalar-taste 
meson is a true Goldstone boson, and for that reason is the valence meson that is usually simulated. 
Here, we calculate the mass of the Goldstone meson in rSZPT. The use of staggered quarks brings 
a new, taste-violating term into the chiral Lagrangian [14]. The two main effects of this poten- 
tial at the one-loop order are the addition of taste breaking corrections to the (staggered) meson 
masses, and the appearance of both a taste-dependent, disconnected vertex and a taste-dependent 
disconnected propagator [14]. 

3.2 Staggered Chiral Perturbation Theory with Electromagnetic Contributions 

Figure 1 shows the 0(p 4 ) contributions to the meson mass. The photon tadpole contributes 
zero in dimensional regularization. The photon sunset is straightforward to calculate, since the 
vertex is taste-conserving. Calculating the contribution from the meson tadpole, Fig. 1 (c), is 
facilitated by using quark- flow diagrams [15]. Figure 2 shows all contributing diagrams considering 
connected and disconnected propagators and vertices. 

The calculation closely follows that of Ref [14], taking care to include electromagnetic effects 
in the LO meson mass. The end result is an electromagnetic contribution from the connected 
diagram, and a non-electromagnetic contribution from the disconnected diagrams that matches that 
of Ref. [14] if in their Eq. (48) is replaced with Bo(m x + m y ). where m x and m y are the masses 

of valence quarks x and y, respectively. 3 




(a) (b) (c) (d) (e) 

Figure 2: Quark flow diagrams corresponding to the tadpole diagram in Fig. 1(c). (a) the diagram with both 
connected propagators and vertices: a "connected" diagram, (b) a disconnected propagator with connected 
vertices and a spectator quark: a "hairpin", (c) a disconnected propagator with connected vertices but no 
spectator quark: a "fancy hairpin", (d) a connected propagator with a disconnected vertex: a "coat hanger", 
(e) both disconnected propagators and vertex: the "coat hanger + bubbles". Note that any of the disconnected 
propagators can include bubble chains. One sees a bubble explicitly in (e), because it is the leading-order, 
disconnected vertex + disconnected propagator diagram. 

3 In Ref. [14], m 2 p+ = #o(m x + m y ). Including electromagnetism, m 2 p+ = 2?o(m x + m y ) + ^Aem- 
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Specifically, for a meson comprised of quark x and anti-quark y, the non-analytic electromag- 
netic contribution is 



S<5 



-2A EM / 1 



31n(M^ 5 /A2)-4 



I6z 2 f 2 V 16 J L % 



(3.5) 



where g X3 , — q x — q y , sea quarks are labeled by a, the sixteen meson tastes are labeled by A x is 
the chiral scale, and £(M 2 ) is the renormalized loop integral 



/ 



d 4 k 



1 



7t 2 k 2 +M 2 



—> £(M 2 ) =M 2 ln(M 2 /A 2 x ) 



(3.6) 



The result in the first line in Eq. (3.5) is from the photon sunset diagram, and that in the second line 
is from the meson tadpole. 

The contributions from Fig. 2 (d) lead to analytic contributions with unknown low-energy 
constants (LECs). We concern ourselves here only with those which are new in the staggered- 
electromagnetic calculation, namely the ones 0(a 2 e 2 ). In principle, one could expect contributions 
from terms like cP^o^o f rom the sea sector 4 , and a 2 (q 2j rq 2 ), and a 2 q 2 y from the valence charges. 
In practice, though, one only has terms °c a 2 {q 2 + q 2 ) and oc a 2 qly. 



4. Utility of Quenched-Photon Simulations 

As was first pointed out in Ref. [13], electromagnetic meson splittings can be extracted at NLO 
in ChPT from quenched-photon simulations with no inherent systematic error due to quenching. 
There are two places where sea-quark charges can appear: in the non-analytic (logarithm) terms 
and in the analytic (LEC) terms. The log terms are completely calculable at NLO. Hence, we can 
put in by hand the difference between the result of a simulation with quenched photons (vanishing 
sea-quark charges) and the value with physical sea-quark charges. For example, the contribution to 
Eq. (3.5) coming from q G ^ can be added back in at the end of the calculation. This is not true 
for the analytic terms, which include, of course, unknown LECs. Nevertheless, the only sea-charge 
dependent LEC terms are 0(me 2 ) and take the form oc [m x + m y ) Y, G tfo- This means they drop out 
of the splittings M 2 ± - M 2 and M 2 ± - M 2 ,. 



5. Conclusions 

We have discussed the origin of electromagnetic effects in quark masses as they arise in a 
typical lattice calculation and shown that the lattice calculation of M^ ± — M^ and M 2 ± — M 2 can 
be used to improve the accuracy and precision of such calculations. For improvement of the calcu- 
lations in Ref. [2], the MILC collaboration has begun staggered-quark, quenched-photon simula- 
tions [3]. In these proceedings, we have described results of the rooted staggered XPT calculation 
needed to extrapolate that lattice data. We have also discussed why quenched-photon simulations 
are sufficient for calculating the needed meson-mass splittings. 

4 The mixed-charge (cross-term) in the sea-sector is zero since (Q) = 0. 
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